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Abstract 

In this paper we consider the possibihty of neutrino mass varying during the evolution of the 
Universe and study its imphcations on leptogenesis. Specifically, we take the minimal seesaw 
model of neutrino masses and introduce a coupling between the right-handed neutrinos and the 
dark energy scalar field, the Quintessence. In our model, the right-handed neutrino masses change 
as the Quintessence scalar evolves. We then examine in detail the parameter space of this model 
allowed by the observed baryon number asymmetry. Our results show that it is possible to lower 
the reheating temperature in this scenario in comparison with the case that the neutrino masses 
are unchanged, which helps solve the gravitino problem. Furthermore, a degenerate neutrino mass 
patten with mi larger than the upper limit given in the minimal leptogenesis scenario is permitted. 
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I. INTRODUCTION 



The baryon number asymmetry of the Universe has been determined to the precision of 
less than 10% after the first year observations of the Wilkinson Microwave Anisotropy Probe 
(WMAP) experiment Q. The reported value of the asymmetry is 

r/i. = - = 6.i;°ixl0-^°, (1) 

where ^ - a„d are the ba,vo„ arrd photon rrumber densities, respective.. 

On the theoretical side, many models have been proposed in the literature [2| to explain 
the small while non-zero dynamically. Among various models, leptogenesis is one of the 
most attractive scenarios 3], where lepton number asymmetry is converted to baryon number 
asymmetry via the(«..).vio>at,n^^haleroninteraet,onsQ . This mechanism has been 



studied extensively in the literature! 



In the minimal scenario of leptogenesis, the standard model is extended by including 
three generations of right-handed (RH) neutrinos, which are Majorana fermions and couple 
to the lepton doublets through Yukawa couplings. The heavy RH neutrinos are produced in 
the early universe by thermal scattering with the primordial thermal bath via the Yukawa 
coupling. In this framework, the cosmological baryon asymmetry is naturally connected 



to the neutrino properties through the seesaw mechanism which is the most natural 
explanation for the tiny neutrino masses observed in the neutrino oscillation experiments 12|. 
As shown in Ref. |3] , the amount of the baryon asymmetry generated in this minimal thermal 
leptogenesis scenario can be characterized by four parameters: the maximal CP asymmetry 
ei, the lightest RH neutrino mass Mi, the effective light neutrino mass m and the quadratic 
mean of the light neutrino masses m, under the assumption that the dominant contribution 
to the lepton asymmetry is given by the decay of the lightest RH neutrino, A^^i. The detailed 
calculations show that to produce the observed baryon asymmetry, it requires Mi > IQ^^GeV 
and Tn % 0.2eV which corresponds to rrii ^ 0.12eV. 

Although it is theoretically elegant and simple, this scenario seems to need a high reheat- 
ing temperature, Tr ~ Mi > O{10^'^GeV), which is only marginally compatible with the 



bound from the gravitino problem, Tr < (10® — 10^^)GeV, in supersymmetric models 
The degenerate solution of the neutrino masses is also strongly disfavored in this scenario. 

n 

Note that if the evidence of neutrinoless double (3 decay with niee = (0.05 — 0.86)eV" 14] (at 
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95% C.L.) is confirmed, a degenerate neutrino spectrum is required. A recent study on the 
cosmological data also showed a preference for degenerate neutrinos with rrii = 0.2eV^^. 
A possible way to accommodate the degenerate neutrinos by the thermal leptogenesis is 
the resonant enhancement of the CP asymmetry when the heavy RH neutrinos are nearly 
degenerate, |M2,3 - Mi| = 0(1^)0,0. 

In the standard ACDM cosmology, the dark energy which drives the universe accelerating 
weights about 73% of the total energy of the universe l|, 3, Il2|- However, the nature of 
the dark energy (DE) remains mysterious. It could be simply a remnant small cosmological 
constant. However, many physicists are attracted by the dynamical solution with a scalar 
field (or multi-scalar fields) like Quintessence. Being a dynamical component, the scalar 
field dark energy is expected to interact with the ordinary matter. In the literature there 
have been a lot of studies on the possible couplings of Quintessence to baryon, dark matter 
and photons 3]- Recent data on the possible variation of the electromagnetic fine structure 
constant reported in 19] has triggered interests in studies related to the interactions between 
Quintessence and the matter fields. 

Recently we have studied various possible interactions, with or without supersymmetry, 



between Quintessence and the matter field of the electroweak standard model 



previous paper 



In a 



2l| we considered a possibility of neutrinos coupling with Quintessence and 



studied its implications on cosmology. In this paper we propose a scenario that the right- 
handed neutrinos in the minimal seesaw model couple to Quintessence and examine its 
implications on the thermal leptogenesis. Other possible irnplications of varying neutrino 
masses in astrophysics and cosmology are discussed in Ref. |2^. With a detailed numerical 
calculation we will show that in this scenario the reheating temperature can be as low as 
about lO^GeV and the degenerate neutrino mass patten with rrii ~ O{0.2eV) will be allowed. 

This paper is organized as follows. In Sec II, we will present the framework and formu- 
lation for the calculation of leptogenesis in our model. We will explain how the interaction 
between the right-handed neutrinos and Quintessence affect rjB iti terms of the four parame- 
ters ei. Ml, m and m. In section HI, we will present the numerical results. In section IV, we 
consider a Quintessence model and propose a specific example of Quintessence to RH neu- 
trinos, then study the interplay between the leptogenesis and the evolution of Quintessence. 
Sec. V is the summary and conclusion of this paper. 
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II. INTERACTING DARK ENERGY SCALAR FIELD WITH NEUTRINOS 



A. The model 



We extend the minimal seesaw model by introducing a new coupling between Quintessence 
scalar field, Q, and the RH neutrinos. The Lagrangian relevant to leptogenesis is given by 

= Ciep + Cq , (2) 

with 

- Ciep = YijLiHNj + ^M,{Q)N^CNi + h.c. , (3) 

and 

Cq = ^d.Qd'^Q - V{Q) , (4) 

where Yij is the Yukawa coupling of the RH neutrinos, Mi{Q) is the Majorana masses of 
the RH neutrinos which is now a function of the value of Q, and V{Q) is the potential of 
the Quintessence. We have taken the basis that the RH Majorana mass matrix is diagonal. 
The form of V{Q) for a specific model and an explicit form of M{Q) will be given in Sec. 
IV. The Quintessence dependent masses of the right-handed neutrinos give rise to different 
Majorana masses at the early epoch from that at the present epoch, which is the key for the 
difference of our model from the usual leptogenesis scenario. 



B. The CP asymmetry 

Generally, the Yukawa coupling Y contains CP violating phases and lead to different 
branching ratios for decays into lepton and antilepton. The asymmetry for A^i decays 



r(iVi ^ LH) - T{Ni LH) 
~ T{Ni LH) + r(A^i ^ LH) 

" -li(ytk7.£'"[(^"^')-rt • M.«M.M,. (5) 

Since the light neutrino masses are generated by the seesaw mechanism and con- 
strained by the neutrino oscillation experiments, the CP asymmetry is also constrained by 



the neutrino data. In the case of hierarchical neutrinos, i.e., JAm^f-j^ ~ O.OSel^ S> 
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^ni^oi ~ O.OOSel^ ^ mi, an approximate upper bound on the CP asymmetry [2l] is 







3 Mn/Am^ 



I I ^ u y atm 

^''^-16^ ' 

with V = llAGeV being the vacuum expectation value of the neutral Higgs boson. For the 



degenerate case, i.e., mi ~ m2 ~ ni-^ ^ m/-\/3 ^ J Arriatrm the upper bound is given by 



^1 ~ 2= — ' V'J 



where m = y mf + m2 + is the quadratic mean of the neutrino masses. One can see that 
the maximal ei is proportional to Mi while inversely proportional to m. In the following 
we always take the upper bound of ei to study the maximal value of tib on the parameter 
space. 

The final lepton asymmetry is determined in general by 

Yl = = K— , (8) 

s 

where s is the entropy density, (?* represents the number of the relativistic degrees of freedom 
at the time when A'^i decays and k < 1 represents the washout effects for the lepton number 
asymmetry in the thermal bath. The lepton asymmetry in Eq. (jHl) will be converted partly 
into baryon asymmetry by the electroweak sphaleron process taking into account the gauge, 



Yukawa interaction and the QCD sphaleron effects 
28/79|lfl, 24|. 



23|, Yb = oYb-l ~ —oyii with a 



C. The washout effect 



The factor k in Eq. (jHI) can be calculated by solving the Boltzmann equations 
At the leading order of ei one has the Boltzmann equations for Y^r^ = rij^^/s and Yb-l 



with 



sHdY^, ^ _ Yr^ 

z dz \ 



- 1 (7fl + + ^lH,t) , (9) 



■ N- 



— ^ = -ei7D [ye,-!)- , (10) 



1 Y 

Iw = 7^7d + 277V + 2jN,t + '2lH,t + tf^7h,s , (H) 
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where 7 is the space time density of the scatterings in thermal equihbrium. In the equations 
above, •jd refers to the reaction density for the process of Ni decays (and inverse decays), 
A^^i ^ LH{LH), 'yH,s and 'yH,t for the processes of AL = 1 scatterings via exchanging Higgs 
boson, NiL{L) ^ t{t)Q{Q) in s channel and Nit{t) ^ L{L)Q{Q) in t channel, 'j^ and 'y^^t 
for the processes of AL = 2 scatterings via exchanging the RH neutrinos, LH ^ LH in 
both s and t channels and LL HH^ LL ^ HH in t channel. 

In Eq. Q, the minus sign on the right-hand side means that all the reactions drive the 
Yatj to its equilibrium value Y^_^. In Eq. flTUI) . the first term on the right-hand side is the 
source term, which gives rise to non-zero Yb-l, while the second term represents washout 
effects that lead Yb-l to zero. 

The consists two parts, the resonant part, '~iN,res and the non-resonant part, '~iN,noni 
i.e., 7Ar = 77v,re.'i + lN,non- The 7d, 7_h-,s, lH,t and '~iN,res are all proportional to the effective 
neutrino mass 8] 

^ {Y^Y),,v' 

while 'jN,non and 'yN,t are proportional to Mim^. These properties are essential in understand- 
ing^he behavor of the solutions the Bolt.maun equations. When S-,„ = + oc 

Miw?' is negligible, i.e., at Mi(m/0.1eV^)^ < 10^^Gel^j8|, the washout factor k is mainly 
determined by m and nearly independent of Mi and m. The effect of m has two different 
aspects. When it is very small (^ 10^'^el^), the RH neutrinos will not be able to reach ther- 
mal equilibrium via the scattering with the thermal bath, so k will be small. Conversely, if 
m is too large the RH neutrinos can be brought into thermal equilibrium rapidly. However, 
the washout effect will also be large, consequently we have a small n. Around m ~ 10~^eV^, 
K reaches its maximal value. When is sizable k will be suppressed exponentially, since 
the washout effect is given by 



with lS.{z) = {Yn^/Y'^^ — l)z/{sH). In this case, k becomes nearly independent of m. 

In short, ei (as well as Y^) increases with^Mi while decreases with m linearly when 
5'~^w is negligible, Mi(m/0.1eV")^ < 10^'^G'el^[8]. However, when 5'~^w oc Mim^ is sizable, 
Mi{m/Q.leVY > lO^^GeV, the final asymmetry Yl will be suppressed exponentially ^. The 
K reaches maximal at m ~ lO^^e^ when the effect of Sjw is negligible. It is shown in Ref. 
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^ that to generate the required value of the baryon number asymmetry, r]B, one needs 
Ml > W^GeV and m < 0.2eV in general. 



D. The effect of Quintessence 

With the interaction between the Quintessence and the RH neutrinos, Mi{Q) is now a 
function of the Quintessence field. Mi at the epoch when it decays will differ from that at 
the present epoch. We introduce a parameter K 

m^) = ^, (14) 

where = Mi{Q^) and Mf = Mi(Q^) denote the Maj orana masses of the RH neutrinos 
at present and at the decay time. For simplicity we take K, as well as Mi, as a constant when 
solving the Boltzmann equations. In section IV, we will show that Q is quite flat during the 
period of leptogenesis for the model we consider and this is a reasonable assumption. We 
then have 

= -Y^Y^v'^ = ml-K , = rff ■ K , = m° ■ , (15) 

where the masses with the superscripts "D" and "0" correspond to the values evaluated at 
the leptogenesis and at the present time respectively. 



III. NUMERICAL RESULTS 



In this section we present our numerical results by solving the Boltzmann equations. 
Besides Mi, m and m we now introduce a new parameter K. We will show that varying 
K makes it possible for Mi to be as low as about IQ^GeV and m can be as large as OAeV 
which corresponds to rrii ^ 

In Fig. ^ we plot the baryon number asymmetry 77^ as a function of K for the hierarchical 
neutrino mass spectrum, i.e., = 0.05e\^. On the left panel, we take = lO^GeV 
for m° = 10~^, 10~^, 10~^eV respectively. One can see from the figure that 7]b increases 

This is the upper bound set by WMAP. Recent SDSs[l3 lowers this hmit to 0.2eV. When neutrino 
masses vary during the evolution of the Universe, however, these cosmological bounds on the absolute 
neutrino masses may be different which is worthwhile studying further. 
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FIG. 1: The baryon number asymmetry rj^ as a function of K for the hierarchical neu- 
trino mass spectrum, i.e., m = OmeV. On the left panel, we take Mf = W^GeV and 
m*^ = 10~^, 10^^, 10~^ey respectively. On the right panel, we take rfi^ = 6 x and 
Aff = 10^, 10^, 10^'^ Gey respectively. The horizontal lines represent tjb = 5.1,6.1,7.1 x 10^^*^ 
respectively. 

with K linearly until ^ 10~^eV, i.e., r]B increases with K until K ^ 10^, 10^, 10^ for 
m° = 10""^, 10~^, lO^^eV respectively. In these regions, the amount of enhancement of ei as 
K gets large dominates over that of the washout effect so that rj^ increases linearly with K. 
However, as the m-^ is about nf ■ K ~ 0.05 ■ lO'^eV^, d-^w becomes sizable then r]B decreases 
exponentially as K. In the region between these two points, is a flat function of K. This 
is due to the effect of K which enhances ei on one hand while decrease n as ^ lO^^el^ 
on the other hand. On the right panel, we fix = 6 x lO^'^eV^, which gives the maximal 
K for the negligible 57vi/|iSj, and take = 10^, 10^, 10^°Ge\^ respectively. One can see 
from this figure that rjB increases with K linearly until (m^/O.leV^)^ ~ lO^^GeV^. The 
different behavior between the two figures is easy to understand if one notices that is 
fixed and does not increase with K on the right panel. 

In Fig. 121 we plot the contour for rjB = 5.5 x 10^^° in the — K plane for the 
hierarchical neutrinos. One can see from this figure that as K gets large, is allowed to 
take a smaller value. For example, at if ~ 30, can be as low as about IQ^GeV . As 
a result of it, the reheating temperature required will be lowered than that in the minimal 
thermal leptogenesis. 

In Fig. El we plot r/^ as a function of K for the degenerate neutrinos, taking wfi = OAeV 
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FIG. 2: The contour oi rjB = 5.5 x 10 ^'^ in the — K plane for the hierarchical neutrino 
spectrum, wfi = 0.05ey. 

{irii ^ 0.23eV) for an example. It should be noticed that m is constrained as mi <m< if 
no fine tunning is takenj^. We take m° ~ ^ 0.23e^ in this case. The figure looks similar 
to that for the hierarchical case. For = (0.23 ■ K)eV ^ lO^^eV, rjB increases linearly 
with K. At the point ijn^ /Q.leVY' ~ lO^^GeV, rjE starts to decrease exponentially as 
K. In the region between, is almost a flat function of K. We notice that for K = 1, 
which corresponds to the case of vanishing Quintessence coupling, no parameter space will 
satisfy the observed rjB- However, for = lO^^GeV, rjB can meet the requirement for 
K ^ 0.5, which is not far from K = 1. One can see that non-zero K makes it possible to 
generate tjb required for the degenerate neutrino mass spectrum. However, in this case we 
have to take K to be smaller than 1 so that n is not suppressed exponentially. This leads 
to a small ei too, which forces us to take a bigger M[^. Therefore, it will not be possible to 
lower and increase Trf simultaneously. 

In Fig. |3]we plot the contour of riB = 5.5 x lO^^'' in the — K plane for the degenerate 
neutrinos. One can see that decreases linearly firstly, then starts to increases slowly. 
At K ^ 0.007, Mf reaches its minimal value, Mf pa 4.5 x lO^^GeV. 
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FIG. 3: The baryon number asymmetry 7]b as a function of K for the degenerate neutrino mass 
spectrum, m° = OAeV. We take = 10^\ lO^^^ XO^^, respectively. The horizontal lines 

represent r]B = 5.1,6.1,7.1 x 10"^'^ respectively. 

IV. QUINTESSENCE MODELS 

In the last section we present our numerical results of leptogenesis taking K as a free 
parameter. In this section we will consider a specific Quintessence model and propose a 
specific form of its coupling to the RH neutrinos. We will study numerically the evolution 
of Q taking into account the back reaction of the neutrino background and calculate the 
factor K. 

In the literature, there have been various proposals for the explicit form of couplings 
in studying the interaction between the Quintessence and the matter field. For example, 
one usually introduces QF^^.F^'^ {Fiiv is the electromagnetic field strength tensor) to study 
the variation of the electromagnetic fine structure constant. However, in an attempt to 
understand the puzzle why the density of dark matter and dark energy are nearly equal 



today, the authors of Ref. [25[ recently consider a model of interacting Quintessence with 
dark matter and in their scenario the mass of dark matter particle depends exponentially 
on the value of Q, m[Q) = ffie'^^^^'^^K 

In this paper, we assume that the coupling between Quintessence and the RH neutrinos 
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FIG. 4: The contour of r]B = 5.5 x 10 in the M/^ — K plane for the degenerate neutrino 
spectrum, mP = OAeV. 



takes a simple form as 

where /? is a 0{1) coefficient. The ratio K is then given as 



(16) 



(17) 



For a numerical calculation of i^', we consider a Quintessence model with the double 

n 

exponential potential |26| 

V = Voie^"^ + e°^) . (18) 

This model has the tracking property for suitable parameters. In the absence of coupling in 
Eq. (fTH|) . the evolution of Quintessence is described by the equations 

SttG 



dQ 



(19) 
(20) 



where H is the Hubble constant, p^, p-y and pg represent the energy densities of matter, 
radiation and Quintessence respectively. We choose the model parameters as A = 100Mp^\ 
a = — lOOMp^^, the initial value of Quintessence field Qi = 1.374Mp; and for the state of 
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FIG. 5: The evolution of Wq and Q as a function of the temperature T for the double exponential 
Quintessence model without the coupling with the RH neutrinos. 



equation, which is defined as 



PQ _ QV2 - V{Q) 



Wq = — = , (21) 

the initial value is Wq- = — 1. We obtain that VLq^^ ^0.72 and the present equation of state 
of Quintessence is PFqq ^ — 1 which are consistent with the observational data. In Fig. El we 
show the evolution of Wq and Q with the temperature T. One can see that Quintessence 
begins to track the background at T ^ IQ'^GeV . 

Taking into account the interaction with the RH neutrinos, we get the equation of evo- 
lution of Quintessence as 



dQ dQ 

The source term in the equation above is given by [2^ 

(3 1 



dViiQ) ^ dM, IM, 



(22) 



(23) 



dQ ^ ~'dQ \ E I MpiTT^' . 
where Ui and E are the number density and energy of the RH neutrinos respectively, () 
indicates thermal average. In the last step of the equation we have taken the Maxwell- 
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FIG. 6: The evolution of Wq and Q as a function of the temperature T for the double exponential 
Quintessence model including the coupling with the RH neutrinos. We take f3 = —2.5 and Mi = 
3.1 X lO^GeV. We take the same definition of Wq as in Eq. (^TJ. 

Boltzmann distribution of the RH neutrinos for simplicity and Ki is the modified Bessel 
function. The energy density of Quintessence is taken the same form as pq = + V{Q) 
in this case. 

We then solve the equation ()22|) numerically, assuming M3 = IOM2 = lOOMi. The 
numerical results are shown in Figs. IHl and [7| In Fig. |Hl we take Mi = 3.1 x 10^ GeV, 
(3 = -2.5, which gives rise to Qo ~ and Qd = l.37AMpi. We then have i^' ^ 31 
and ^ lO^GeV, corresponding to the case we considered in the hierarchical neutrino 
spectrum. In Fig. Owe choose the parameters Mi = 3.15 x lO^^GeV^ and (3 = 3.62. We find 
the values of Qq and Qn are almost the same as the above case. We then have K ^ 0.007 
and Mf ^ 4.5 x corresponding to the case that satisfy our requirement for the 

degenerate neutrinos. 

Comparing Figs. IHlandElwith Fig. one can see that the interaction of Quintessence with 
the RH neutrinos does change the early behavior of the Quintessence field Q and its equation 
of state, however, does not change the tracking properties of this model. Furthermore, the 
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FIG. 7: The evolution of Wq and Q as a function of the temperature T for the double exponential 
Quintessence model including the coupling with the RH neutrinos. We take (3 = 3.62 and Mi = 
3.2 X W^^GeV. We take the same definition of Wq as in Eq. (ETJ). 

value of Quintessence field Q remains a constant in this model until T ~ lO'^GeV which 
satisfies our assumption for a constant K during the period of leptogenesis. 



V. SUMMARY AND CONCLUSION 



The explanation of the baryon number asymmetry of the Universe remains a challenge 
for cosmological and particle physics. Among various theoretical models, leptogenesis is one 
of the most attractive scenarios. Especially, the minimal thermal leptogenesis only needs 
the minimal extension of the standard model, while the extension seems necessary in order 
to explain the neutrino oscillation experiments. 

The minimal thermal leptogenesis is closely connected with the properties of the light 
neutrinos. Detailed studies show that, to generate the observed baryon asymmetry, the 
constraints Mi ^ lO^^GeV^ and rrii ^ O.ley should be satisfied. These results disfavor the 
degenerate neutrino mass spectrum strongly and require a high reheating temperature. 
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We investigate in this paper the imphcations on the minimal thermal leptogenesis of a 
new scenario that the RH neutrinos couple with the Quintessence, whose evolution causes 
the masses of the RH neutrinos vary. We study the interplay of this coupling between 
the evolution of Quintessence and the leptogenesis. By solving the Boltzmann equations 
numerically we show that Mi, at the decay time, can be as low as about lO^GeV^ and mj, 
at the present epoch, can be as large as ~ 0.2eV. 
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